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To implement a set of universal quantum logic gates based on non-Abelian geometric phases, it is a conven¬ 
tional wisdom that quantum systems beyond two levels are required, which is extremely difficult to fulfill for 
superconducting qubits and appears to be a main reason why only single-qubit gates were implemented in a 
recent experiment [A. A. Abdumalikov Jr. et al., Nature (London) 496, 482 (2013)]. Here we propose to realize 
nonadiabatic holonomic quantum computation in decoherence-free subspace on circuit QED, where one can use 
only the two levels in transmon qubits, a usual interaction, and a minimal resource for the decoherence-free sub¬ 
space encoding. In particular, our scheme not only overcomes the difficulties encountered in previous studies, 
but also can still achieve considerably large effective coupling strength, such that high fidelity quantum gates 
can be achieved. Therefore, the present scheme makes realizing robust holonomic quantum computation with 
superconducting circuits very promising. 

PACS numbers: 03.67.Lx, 42.50.Dv, 85.25.Cp 


I. INTRODUCTION 

Under adiabatic cyclical evolution, a quantum system ac¬ 
quires a phase factor, which consists of both dynamical and 
geometric components. When the eigenstates of the sys¬ 
tem are nondegenerate, the geometric component is the well- 
known Berry phase Ol]. For the degenerate case, it is a unitary 
aerator acting on the degenerate subspace, i.e., holonomy 
11^ . As geometric phases are determined by the global prop¬ 
erty of the evolution path, the geometric method of quantum 
computation has been shown to possess some built-in noise- 
resilience features lEHU. In general, the holonomies do not 
commute with each other, and thus can be used to construct 
a universal set of quantum gates lEMl, i.e., the holonomic 
quantum computation (HQC). 

On the other hand, the adiabatic method of quantum com¬ 
putation intrinsically leads to long gate operation time, which 
may be comparable with the lifetime of typical qubits ll^l2^ . 
This motivates research on quantum computation based on the 
nonadiabatic geometric phases. Recently, nonadiabatic HQC 
has been proposed using three-level A systems ll^ with the 
experimental implementation of some elementary gates & 
1^ . However, the excited state is resonantly coupled when 
implementing the quantum gates 1^ . and thus, its limited 
lifetime is a main challenge in practical experiments. Note 
that this limitation may be avoided in experiments in Refs. 

and 0 because they use the three magnetic states of 
a nitrogen-vacancy center in a diamond. However, for a su¬ 
perconducting transmon qubit, this limitation does exist, and 
recent experiment has verified only sin gle- qubit gates 12^ . 
The energy levels of a transmon qubit are in a ladder 
shape, and the anharmonicity is small, which limits the cou¬ 
pling strength between neighboring levels to the order of 10 
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MHz in order to individually address the interactions 
Therefore, even with newly demonstrated good coherent times 
of multilevels in the transmon qubit 1132] . the implementa¬ 
tion of a nontrivial two-qubit holonomic gate, which needs 
much more complicated cavity-induced interaction between 
two three-level systems is still very challenging. Alter¬ 
natively, there are schemes using circuits more complicated 
than transmons to mimic a multilevel system ||3,[lal^- How¬ 
ever, this will inevitably introduce additional noises from the 
environment because more circuits and control elements are 
needed. 

Meanwhile, many efforts have also been made to com¬ 
bine H QC w ith the decoherence-free subspace (DFS) encod¬ 
ing JUflllt]. HQC in DFS Isbl - lddl] can consolidate both the 
noise resilience of the encoding and the operational robust¬ 
ness of holonomies. As for transmon qubits, this protocol 
is much more difficult to implement because it requires at 
least two transmon qubits to encode a logical qubit, and thus, 
more complex interactions among qubits are needed, even in 
the single-qubit-gate case. Moreover, previously proposed 
schemes based on HQC in DFS usually need at least three 
physical qubits to encode a logical qubit; note that the use of 
two physical qubits is a minimum resource needed, as in Ref. 
for geometric entangling gates. 

In this paper, we propose to implement a nonadiabatic HQC 
in DFS with a typical circuit QED setup. Our scheme avoids 
the above-mentioned difficulties. First, only the two energy 
levels of the transmon qubits are involved. Second, for the 
single-qubit case, our implementation relies solely on the ef¬ 
fective resonate qubit-cavity interaction. Meanwhile, for the 
two-qubit case, only a conventional detuned interaction is re¬ 
quired, where the detuning between a transmon qubit and the 
cavity is fixed, and thus, we have plenty freedom to avoid 
the limitation due to the small anharmonicity of transmons. 
Third, we use two transmons to encode a logical qubit, which 
is the minimal resource for the DFS encoding. Therefore, our 
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FIG. 1. (Color online) Illustration of the setup of our scheme, (a) 
The coupled two-transmission-line-resonator system for universal 
holonomic quantum computation in a decoherence-free subspace, in 
which each resonator has two transmon qubits to encode a logical 
qubit, (b) The qubit and cavity coupling configuration. 


Wqj(f) = Wq -F Ej sm{ujjt - If j). (2) 


This modulation may effectively turn the qubit’s sideband on 
resonance with the cavity frequency. This can be clearly seen 
by moving to the rotating frame defined hy U = Ui x U 2 , 


Ui = exp 
U 2 = exp 


—i 


(^cr| + Wco'l'a^ t 


ia^^ cos{ujjt - fj) 


( 3 ) 


with aj = Ejlujj, and the transformed Hamiltonian is 


scheme presents a promising method for HQC on supercon¬ 
ducting circuits. 


II. THE SETUP AND INTERACTION 


= {gUjCL^ + H.c.) U 

= (gcTj-a^-F H.c.) U 2 (4) 

= gajo )exp [iaj cosiuijt — ipj)] + H.c.. 


The setup we consider is illustrated in Fig. 1(a), which 
consists of two coupled transmission line resonators (one¬ 
dimensional cavities). In each cavity, there are two transmon 
qubits ll^ capacitively coupled to it. This coupled system is 
the building block of our scheme, and the two transmon qubits 
encode a logical qubit. For the case of single-logical-qubit 
gates, we consider only the dynamics within a cavity, while 
the coupling between the two cavities is switched on when 
implementing the two-logical-qubit gate. Furthermore, one 
may repeat this block to construct a one-dimensional chain of 
logical qubits; that is all the parameters of the odd and even 
qubits are the same as those of the first and second qubits, 
respectively. 

We first detail our setup for the single-qubit case. A trans¬ 
mon qubit is composed of two identical Josephson junctions 
in a loop configuration and is shunted by a large capacitance. 
The transmon qubit is quantized and its lowest two energy lev¬ 
els can be used to construct our physical qubit states with the 
effective Hamiltonian iTqj = a;qjtT|/2 {h = 1), where Wqj is 
the energy splitting of the transmon qubit and cr^ is the Pauli 
matrix of the jth transmon qubit in its eigenbasis. For ty pica l 
values of experimental parameters, Wqj ~ [4,10] GHz 114^ . 
The transmon qubits are located at the voltage antinodes of 
the relevant cavity mode, = oJca^a, where Wc, a, and 
are the frequency, annihilation, and creation operators of the 
cavity, respectively. The coupled system is described by ll4^ 


Hic = H, 


i=i 


[Hqj + gicTjCL^ + crja) 


( 1 ) 


where g is the qubit-cavity coupling strength (assumed to be 
real), aj is the transmon lower operator, and crj = (aj)^. Here 
we consider the case of A = (wc — Wq) ^ ; that is, the qubit- 
photon interaction acts perturbatively. 

To get a resonate interaction between a selected transmon 
qubit and the cavity, the qubit is biased by an ac magnetic flux, 
which will introduce periodical modulation jidl] of the qubit 


Using the Jacobi-Anger identity of 

exp[iQ;j cos{u!jt — fj)] 

00 

= ^ i'^Jm{aj)eyi-p[im{LOjt- fj)], 
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and J-m{ctj) = with Jm{oij) being Bessel 

functions of the first kind, the transformed Hamiltonian re¬ 
duces to 


Hd,i = gJo{aj){ajo)e'-^^ + H.c.) (5) 

00 

-F ga^a) ^ + H.c.. 

m—1 

When ujj = A, the effective resonate qubit-cavity interac¬ 
tion will be in the form of 


Hij = gjCTja'^ + H.C., (6) 

where gj = gJi{aj)exp{iipj + 7r/2), we have applied 
the rotating-wave approximation by neglecting the oscillating 
terms, and the smallest oscillating frequency is A, i.e., 


iTosc = guja^ 


H.C.. 

( 7 ) 


In this way, we can have full control of the coupling strength 
gj by varying the externally driven ac magnetic flux, i.e., by 
controlling the amplitude Sj and phase ipj. Meanwhile, as¬ 
suming that the anharmonicity of the transmon is the same as 
A, as the transition frequency of |1) o |2) is ojq — A, the 
third level of the transmon can only couple dispersively with 
the detuning 2A, as shown in Fig. 1(b). Finally, note that 
the resonate interaction in Eq. (|6l), for the case with more 
than two qubits, has the conserved quantity of total excitation 
N = X]j=i with ric being the photon number in 

the cavity. 
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III. SINGLE QUBIT GATES 


We now proceed to deal with the holonomies for single¬ 
qubit gates in DFS. Hereafter, to avoid confusion, we refer to 
our physical transmon qubits as transmons and logical qubits 
as qubits for short. As the transmons are placed in the same 
cavity, they can be treated as interacting with the same cavity- 
induced dephasing environment. The DFS we consider here 
is the subspace of 

= {|100), 1010), 1001)} = {|0)i, |l)i, \E)l}. (8) 

where the subscript L denotes the states belonging to the log¬ 
ical qubit and | 100 ) = |l)i ® | 0)2 ® | 0 )c, i.e., they denote the 
states of the first and second transmons, and the cavity, respec¬ 
tively. Also, we use the cavity as an ancillary, and thus, only 
two transmons are needed to encode a logical qubit. Note that 
the DFS is identical to the subspace of iV = 1, which ensures 
that the quantum dynamics will not go out of subspace Si. 

In this encoding, the Hamiltonian of the quantum system 
consisting of two transmons, i.e., j € { 1 , 2 }, resonantly cou¬ 
pled to a cavity reduces to 

Hi=gi\E)L{Q\+g2\E)L{l\+^.c. 

= a (^sin^e*‘"|F;)L(0|-cos^|F;)L(l|+H.c.^ ,(9) 

where a = Ji{ctiY + 'h{ci 2 )^ is the effective Rabi fre¬ 
quency, tan( 0 / 2 ) = Ji{ai)/Ji{a 2 ), and (/? = <^i — — tt. 

In this case, we construct a A-type Hamiltonian in the DFS 
with only resonate transmon-cavity interaction, from which 
an arbitrary single-qubit holonomic gate can be obtained. It is 
worth pointing out that the two transmons do not have to pos¬ 
sess the same frequency, as we may use two externally driven 
fields with different frequencies to bring them in resonance 
with the cavity. 

In the dressed-state representation, the Hamiltonian in Eq. 
@ can be viewed as indicating that state \E)l couples with 
only the ’’bright” state \b) = sin |e“®‘^| 0 ) — cos || 1 ), while it 
decouples from the ’’dark” state \d) = cos || 0 ) -f sin |e®‘^|l). 
Under the action of Hi, the dark and bright states evolve ac¬ 
cording to 

\Mt))L = Ui{t)\d) = \d), 

iV^ 2 (f))L = Ui{t)\b) = cos(^if)|&) - isin(^if)|£')i.^ 

When the condition = tt is satisfied, the dressed 

states undergo a cyclic evolution as |' 0 i(ri))(^/)i(Ti)| = 

\ipi{0)){ipi{0)\. Under this condition, the time evolution op¬ 
eration on the subspace (Id), |&), li?)^} is given by 


Ui{ti)= ^ 


.JV'i(0))(r/^j(0)l, (11) 




where T is the time-ordering operator and Aij{t) = 
In particular, when the condition iTij(f) = 
{ipi{t)\Hi\'ipj{t)) = 0 is met, which means that there is no 
transition between the two time-dependent states, the evo¬ 
lution satisfies the parallel-transport condition. Therefore, 
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FIG. 2. (Color online) The performance of the proposed Hadamard 
gate. (a) Qubit-state population and fidelity dynamics of the 
Hadamard gate as a function of dimensionless time (b) Max¬ 

imum fidelity for different cavity decay rates k (in units of k) and 
A/p = 10, with p/(27r) being 50 and 100 MHz for the blue circles 
and red triangles, respectively. 


the geometric nature of the operation is originated from the 
structure of the Hamiltonian instead of the slow evolution in 
the adiabatic case. Under these two conditions, in the logi¬ 
cal qubit subspace {|0)l, |1)l}, the nonadiabatic holonomic 
gates can be realized as 


U{e,g,) 


COS0 sin0e \ 
sin 06*'^ — cos 0 J ’ 


( 12 ) 


where 0 and tp can be tuned by choosing approximate param¬ 
eters of the externally driven ac magnetic fluxes. Therefore, 
an arbitrary single-qubit gate can be achieved. 

The performance of the gates can be evaluated by consid¬ 
ering the influence of dissipation using the quantum master 
equation; 

p = i[p, Hi + H'l] + |-C(a) 

-\-—C{ai+a2) + —C{a\+a\), (13) 

where p is the density matrix of the considered system; 
C{A) = 2ApA^ — A"^Ap — pA^A is the Lindblad operator; 
and K, Fi and F 2 are the decay rate of the cavity and the de¬ 
cay and dephasing rates of the qubits, respectively. We have 
assumed that the decay and the dephasing rates of the two 
transmons are the same. We consider the Hadamard gate to 
be a typical example, where 0 = 7 r /4 and p = Q. To make 
the total coupling strong, we may choose Ji(ai) ~ 0.207 
and Ji(a 2 ) = 0.5, which corresponds to Ji(q;i)/Ji(q; 2 ) = 
0.414 ~ tan(0/2). This can be achieved by modulating 
cti = ei/A ~ 0.4236 and a 2 = £ 2 /A ~ 1.2068; note that 
tuning the coupling strength in such a way has been experi¬ 
mentally demonstrated 114^ . We may choose A = 27r x 500 
MHz, 5 = 27r X 50 MHz, and thus ~ 0.54p = 27r x 27 
MHz. To verify the approximation when obtaining Eq. (|6l), 
we also include some oscillating terms in the simulation. As 
Jo(Qfi)/J 2 (Q!i) > 40, the J 2 {oii) term can be safely ne¬ 
glected for transmon 1. Meanwhile, for transmon 2, this ratio 
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is also larger than 4. Therefore, for demonstration purposes, 
we include only the Jq (aj ) terms in Eq. O for the two trans- 
mons in our numerical simulation, i.e., in Eq. (fOl) . 

2 

H[ = + H.c.] . (14) 

i=i 

Quality factors of about 10® have been reported for a cavity 
with frequencies ranging from 4 to 8 GHz 1^ . and thus, the 
cavity decay rate k is on the order of kilohertz. Relaxation 
and coherence times of 44 and 20 /rs are reported for a 
planar transmon, which corresponds to Ei ~ 27r x 8 kHz and 
r 2 ~ 27r X 3.5 kHz. As k, Ei and E 2 are all on the same 
order of magnitude, for simplicity, we treat them as if they 
were identical and set Ei = E 2 = k = 27r x 10 kHz. Sup¬ 
pose the qubit is initially in the state |0 )l; we evaluate this 
gate by the qubit-state population and the fidelity defined by 
Fi = {ipflplipf), with 1'!^/) = (|0) -f |1 ))l/V 2 being the 
ideally final state under the Hadamard gate. We solve nu¬ 
merically the master equation ( fTST l with H[ being given in 
Eq. (|7]); as shown in Eig. |2ja), we obtain a very high fidelity 
Fi ~ 99.8% at t = tt/^i ~ 18.5 ns. Meanwhile, the higher 
energy levels will also be involved during the evolution; we 
next consider the effect of the third level of the transmons, 
which is the one closest to Wq. Note that more higher ex¬ 
cited levels will not directly influence the transmon qubit sub¬ 
space. Eor this level, we assume that the anharmonicity of the 
transmon is the same as A, and thus, the cavity-induced cou¬ 
pling between the second and third levels when obtaining the 
Hamiltonian in Eq. (|6l) is dispersive, i.e., oscillating with fre¬ 
quency 2A. Therefore, its contribution to the infidelity of the 
quantum gates is similar to that of the Hamiltonian in Eq. (|7]). 
Assuming that the coupling strength is V^g, the infidelity that 
results from this third level is 0.06% from our numerical sim¬ 
ulation. In addition, when the transmons are incorporated into 
the cavity, its decay rate will be enhanced. Therefore, we also 
investigate the influence of the increase of the cavity decay on 
the gate fidelity; as shown by the blue circles in Eig. |2tb), for 
k' = 30k, we can still get Fi ~ 95%. Meanwhile, the larger 
detuning implies g can be bigger, and thus, the fidelity can be 
higher. As shown by the red triangles in Eig. |2b), we verify 
this by choosing the detuning and transmon-cavity coupling 
strength to be 2g and 2A, respectively. 


IV. TWO-QUBIT GATE 

At this stage, we turn to the implementation of a nontrivial 
two-qubit gate. To avoid the cross talk between transmons, 
we consider the scenario of two coupled cavities, with each 
one having two transmons to encode the logical qubit. In this 
case, a six-dimensional DES exists 

5*2 = {|00)l = 1100100), |01)l = llOOOlO), 

|10)l = 1010100), |11)l = 1010010), 

|ai) = 1110000), |a2) = |000110)}, (15) 

where |10100) = |l)i| 0 ) 2 | 0 )ci|l) 3 | 0 ) 4 | 0 )c 2 ; |ai) and |a 2 ) are 
two ancillary states, both of which have two excitations within 


a logical qubit and thus will not be affected in the single-qubit 
cases (in the N = 1 subspace). 

To obtain a nontrivial two-qubit gate, we need to induce the 
interaction within two transmon pairs, i.e., transmons 2 and 3 
and transmons 2 and 4, but to avoid the interaction between 
other pairs. To achieve this, we consider that the interaction is 
induced by exchanging virtual photons between the two cou¬ 
pled cavities, the coupling of which is 

iTcc = X{a\a 2 + aial), (16) 

and the frequencies of the delocalized held modes Pi = 
(oi — a 2 )/\^ and P 2 = (oi -f a 2 )lV 2 are shifted lld^ from 
the bare cavity frequency as wi = Wc — A and UJ 2 = uJc + A, re¬ 
spectively. In addition, to avoid cross talk between the two in¬ 
teracting pairs, they have different detunings. The setup is de¬ 
tailed as follows. Eor transmons 3 and 4, which are located in 
the second cavity, we do not apply externally driven on them 
and set S = (wq^ — Wc) = (wc — <^^, 4 ) = 27r x 150 MHz, 
(73 = —gi = <72 = <? = 27r X 30 MHz, and A = 25. As for the 
third energy level, the anharmonicity of the transmon is ap¬ 
proximately A = 27r X 500 MHz, and thus, the higher levels 
can be safely neglected. Then, in the bare cavity frequency, 
the interaction Hamiltonian reads 

= ga\{a^e-^^* - + H.c.. (17) 


Meanwhile, for transmon 2, we modulate the frequency of the 
driven ac magnetic flux as = 25. In the rotating frame, the 
interaction Hamiltonian reduces to 


Hm 2 = ga\a2 fjo(/3)e*‘'* + 


H.c.,(18) 


where (3 = 4’ = ‘f 2 ~ ^/2- 

However, in the presence of Hcc, the frequencies of the 
delocalized cavity held modes Pi = (ai — a 2 )l'j 2 and 
P 2 = (oi + a 2 )l \/2 will be renormalized from the bare cav¬ 
ity frequency as wi = — X and a ;2 = Wc + A, respectively. 

Then, in the interaction picture, the held mode operators will 
be renormalized as 

4 ^ . (19) 

Therefore, the total interaction Hamiltonian reads 


Hint — H'lntX 

^ 9 


Pint2 

4e*(^+^)* -f 


H.C., (20) 


where h\ = p|[Jo(/3)(T2 — < 74 ] + Pi[Ji(,8)e *'^4 “ 4] ™tl 

h\ = Pj [Ji(/ 3 )e®‘^(T 2 + 0 - 3 ] -I- Pi[Jo(/ 3)4 + 4]- Assuming 

that {A —15, 25} ^ glV^, the above interaction can be treated 
as if it has two independent interaction channels that oscillate 
with distinctly different frequencies, with the cross talk be¬ 
tween them being suppressed by a frequency difference of 25. 
When 25 ^ gjV^, the cross talk can be safely neglected. The 
effective Hamiltonian of the total interaction is 


Heff = 77 


Ji(/3)e*‘^cr24 “ 


H.C., (21) 
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FIG. 3. (Color online) Qubit-state population and fidelity dynamics 
of the (7(7r/4, 0) gate as a function of dimensionless time ^2f/7r. 


where rj = — <5^) and the Stark shift term has been 

neglected. In order to turn off this coupling we may 

modulate the coupling strength to be time dependent as X(t) = 
2 Acoswf, as recently demonstrated experimentally 0 ^ 
and the two cavities have a frequency difference of a; = | Wci — 
U!c 2 \, with ujci and ojci being the resonant frequencies of the 
first and second cavities, respectively. 

In subspace S 2 , the Hamiltonian in Eq. (I 2 TI 1 becomes 




sin^e*‘^(|ai)i(00| -F |ll)L(a2|) 




-cos-(|ai)L( 01 | -F 110)^(02!) +H.C. 


, ( 22 ) 


where tan(i?/2) = Ji(/3)/Jo(/3) and the effective Rabi fre¬ 
quency ^2 = vVJ q{P)^ + The effective Hamilto¬ 

nian in Eq. (|2^ can be divided into two commuting parts as 
H 2 = Ha + Hb, with 


Ha = ^2 


sin-e*'P|ai)L 


cos-|ai)L( 01 | 


H.C., 


Hb=^2 


sin^e *'^|a2)c(ll| - cos^|a2)c(10| 


H.c. 


When ^ 2 T 2 = tt, the evolution operator in our logical qubit 
subspace reduces to 




/ cost? 
sin t?e*'^ 


0 

V 0 


sint?e-*'^ 0 

— cos t? 0 

0 — cos t? 

0 sin t?e*'^ 


0 

sini?e“®'^ 
cost? J 

(23) 


We can see that the gate in subspace {|00), |01)} is different 
from the one in subspace {|10), |11)}. Therefore, in general, 
this is a nontrivial two-qubit gate. 

Eor example, when t? = 7 r /4 and </> = 0, it reduces to 



/I 

1 

0 


( 4 ’^^) 72 

1 

-1 

0 

0 

0 

0 

-1 

1 


\0 

0 

1 

1 / 


In this case Ji{P)/Jo{0) = tan( 7 r/ 8 ), which leads to /3 ~ 
0.77, Jo(/3) = 0.86, Ji(/3) = 0.36, and thus ^2 ^ x 6.2 
MHz. In addition, Ji(0.77)/J2(0.77) > 5 and J2(0.77) > 
J„(0.77) for n > 2. Moreover, the smallest oscillating fre¬ 
quency is 3(5 when n = 2, and thus, the coupling between 
two transmons induced by the J 2 ( 0 . 77 ) term is much smaller 
in strength compared with that induced by the Ji(0.77) term. 
Therefore, the higher-order terms are neglected in Eq. ( fTST l. 
Eor the initial state |01), we simulated the performance of 
this gate using the master equation with the total Hamiltonian 
Hi = i7cc + ?Tinti +?Tint 2 inEqs. (fTfil l. (fTTl i. and (fTSl l. as shown 
in Eig. 3, where a high fidelity of F 2 ~ 94.5% can be reached 
with the decay rates being the same as in the single-qubit case. 
The infidelity mainly comes from both the decoherence of the 
system and the validity of the effective Hamiltonian in Eq. 
(EB, which are 3% and 2.5%, respectively. Meanwhile, con¬ 
sidering the influence of the third level will additional intro¬ 
duce about 0.5% infidelity. In addition, when the anharmonic- 
ity is smaller, e.g., 27r x 300 MHZ, the two-qubit fidelity will 
be 93.1% with p = (5/5 = 27r x 20 MHz. Einally, as transmon 
3 has a different frequency than the others, its coupling should 
not be the same as the others. This problem can be solved by 
redefining tan('!?/2) = 53 Ji(/3)/[pJo(^)]. 


V. CONCLUSION 


In conclusion, we have proposed to implement HQC in 
DES with the typical circuit QED, in which only two levels 
from conventional transmon qubits are required. 
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